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Based on first-principles calculations, we demonstrate that magnetism impedes the formation of 
long chains in break junctions. We find a distinct softening of the binding energy of atomic chains 
due to the creation of magnetic moments that crucially reduces the probability of successful chain 
formation. Thereby, we are able to explain the long standing puzzle why most of the transition- 
metals do not assemble as long chains in break junctions and provide thus an indirect evidence that 
in general suspended atomic chains in transition-metal break junctions are magnetic. 



One-dimensional systems, realized experimentally as 
suspended monatomic chains in break junctions (BJs), 
have altered our conceptional view on atomic scale junc- 
tions. For example, due their enhanced tendency to 
magnetism they bear high potential in the field of spin- 
tronics by combining the possibility to probe, control 
[H, [H, and switch the magnetic state by spin-polarized 
electrical currents. Recent first-principles calculations 
support these expectations and report on sizeable mag- 
netic moments [3|, 0, S H J3 an d giant magneto- 
crystalline anisotropy energies 0j 0, 0] m suspended and 
free-standing transition- metal (TM) monowires (MWs). 
While the formation of long atomic chains of selected 
TMs and conductance quantization have been experi- 
mentally demonstrated [g, @, > any conclusive evidence 
of magnetism in chains is still missing. 

This lack of evidence seems even more surprising as 
measurements on chains deposited on surfaces show uni- 
vocal signatures of local magnetic moment 111, LL2J . Mea- 
suring the magnetoresistancc in BJs would serve as a 
proof for magnetism in atomic-sized contacts. In ma- 
terials, where long MWs can be successfully suspended 
as reported for Ir and Pt BE El, the leads are non- 
magnetic, which prevents the pinning tip-magnetization 
and thus the analysis via magnetorcsistivity measure- 
ments. On the other hand, so far it has not been shown 
that breaking contacts of magnetic 3d-TMs results in 
one-dimensional structures beyond point contacts. 

An alternative approach to prove that monatomic 
chains are magnetic is the search for half-integer con- 
ductance originating from 100% spin-polarized conduc- 
tance channels. Although Rodrigues et al. [l5| reported 
on such a half-integer conductance for Co, Ni and Pt BJs, 
it was shown both experimentally [l6| and theoretically 
17J that not only magnetic ordering but e.g. also the ad- 
sorption of H2 can lead to similar values of conductance. 
In other words, the presence of half-integer conductance 
is not a unique attribute of full spin-polarization and thus 
cannot serve as a proof of it. 

Because of these uncertainties, in this letter we take 
a new path to address the emergence of magnetism in 



BJs. We base our study on the most fundamental and 
easiest accessible experimental quantity: the probabil- 
ity for successful chain formation of a given material it- 
self. The trend arising from numerous BJ experiments 
is, that monatomic chain formation is most probable for 
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late 5c? TMs as well as Ag and Au 
analyze the role of magnetism for the chain formation 
in BJs of 3d, Ad and hd TMs we apply a recently de- 
veloped material-specific theoretical model [l9[ for the 
formation of long monatomic chains which operates in 
terms of parameters extracted from ab initio calculations. 
By explicitly including and excluding magnetic exchange 
interactions, we prove that magnetism significantly sup- 
presses the chain formation due to a substantial reduction 
of the chain hardness expressed in terms of the maxi- 
mally sustainable break force. Comparing our results to 
experimental findings, we are able to provide an indirect 
evidence that chains in BJs arc indeed magnetic. 

We will briefly recall the model which allows to inves- 
tigate the formation probability of suspended monatomic 
chains under tension in BJs [191 ] before applying it to both 
magnetic and non-magnetic chains. In this model, chain 
formation succeeds if the criteria for stability and pro- 
ducibility are met. The criterion for stability addresses 
the rupture of the chain via breaking of the bond be- 
tween two neighboring chain atoms. The criterion for 
producibility is concerned with the chain elongation that 
is composed of two processes: At first one atom has 
to be extracted out of the lead into the chain. This 
transfer leads to a reduction of the coordination of this 
particular atom and consequently additional external en- 
ergy is required, which we account for by the difference 
A-ELead = %w — -^Lead of the energy of cohesion for an 
atom in the lead, -B^ead, and in the chain, Emw, both at 
equilibrium distance. This energy can be stored mechan- 
ically and is released in the second process of relaxing all 
chain's bonds to a smaller inter-atomic distance after an 
additional atom has joined the chain. These competing 
contributions to the total energy of the system determine 
whether it is energetically favorable to grow the chain by 
one atom or not. In order to apply both criteria only 
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FIG. 1: (color online) Calculated MW energy £(d) (circles) as 
a function of inter-atomic distance d (throughout the paper 
given in a.u.=0.0529nm) for both non-magnetic (NM) (black) 
and magnetic (M) state (green) are shown for (a) Fe and 
(b) W. The Morse-fit (lines) provides a universal fit to these 
points well below the required accuracy of about 100 meV for 
chemical bonding and can be characterized by the equilibrium 
inter-atomic distance do and the inflection point d. £(d) and 
the magnetic moments M in (c) and (d) correspond to the 
FM and AFM ground state for Fe and W, respectively. 



two quantities have to be known: The binding energy 
£(d) of a MW atom as a function of inter-atomic dis- 
tance d and the cohesive energy difference Ai?Lcad- We 
will show that both parameters depend on the formation 
of magnetic moments. 

In order to determine these parameters we carried 
out spin-polarized and non-spin-polarized ab initio cal- 
culations in the generalized gradient approximation 
(GGA) [2(| to the density functional theory for selected 
3d, Ad, and 5d TMs, employing the full-potential lin- 
earized augmented plane-wave method (FLAPW) for 
one-dimensional (ID) systems plj . as implemented in the 
FLEUR code [22[ . We used the bulk version of this code to 
calculate the cohesion energy [23| . Here we considered in 
all cases the true magnetic (3d) or non-magnetic (4d, 5d) 
bulk ground state as the reference configuration. For cal- 
culations of the bare MWs in all cases we considered the 
non- (NM), ferro- (FM) and antiferromagnetic (AFM) 
order and included basis functions with plane waves up 
to fcmax = 4.4 a.u." 1 and used 64 k points in one half 
of the ID Brillouin zone. We calculated all 3d MWs in 
the scalar-relativistic approach, while spin-orbit coupling 
was added for all 4d's and 5d's. 

The calculated spin moments M(d) and the binding 
energy £(d) are shown in Fig. 1 for Fe and W MWs as 
examples. Both chains exhibit sizeable magnetic mo- 
ments and the binding energies for the non-magnetic 
(NM) and magnetic (M) case differ substantially. This 
is true in particular for large distances, where the en- 
ergy difference £ sp (d) = £nm(gO — £iA{d) approaches the 



spin-polarization energy £ sp (oo) of an isolated atom. A 
closer look reveals, that this overall tendency to mag- 
netism in ID chains results in finite magnetic moments 
even in 5d-TMs already at the NM equilibrium distance 
dp M or upon small stretching leading to a magnetic ex- 
pansion of e.g. d^ - 4™ = 0.2 a.u. for W (Fig. 1(a)) 
and being more pronounced e.g. for Fe (Fig. 1(b)), in 
agreement with Ref. 0, [2f|. Overall, the binding po- 
tential energies of M and NM chains differ not only by 
spin-polarization energy and a constant shift in do, but 
also their slopes are crucially different (Fig. l(a),(b)). 

In order to analyze the binding energy quantitatively, 
we fit a Morse-potential 



£(d) = £(oo) • (l 



_ e -l(d-d ) 
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to the discrete set of calculated energy points 
(Fig. l(a),(b)) for the different magnetic states. Besides 
the equilibrium distance do, the Morse potential can be 
characterized by the following two physically transparent 
parameters: The inflection point d = do + ln2 /j and the 
break force Fq = 7 • £{oo)/2, which is the maximal slope 
F(d) of the potential. Together with A£"Lead> which we 
evaluate for a given close-packed surface (Fig. 2(a)), these 
four quantities constitute a minimal basis for a realistic 
description of the chain formation process, providing us 
with an accurate and continuous representation of the 
binding energy curve. 

As we see from Fig. 2, magnetism has a very differ- 
ent impact on these four quantities: The cohesion energy 
difference AE'Lcad (Fig. 2(a)) of 4d and 5d TMs remains 
almost unchanged for magnetic and non-magnetic chains, 
as these elements are non-magnetic at surfaces and de- 
velop only small magnetic moments at do- For 3d TMs 
the situation is different: in the wire the magnetic mo- 
ments at do are already close to the saturated values and 
larger than for atoms at surfaces, which leads to a size- 
able gain in spin-polarization energy and thereby to a 
considerable reduction of A£"Lead |26j. The position of 
d with respect to do (Fig. 2(b)), a quantity crucial for 
the stability of the wires, reveals only minor changes of 
less than 5% between calculations with and without spin- 
polarization. The crucial impact of magnetism is the soft- 
ening of the binding energy curve, which leads to a sig- 
nificant reduction of the break force Fo for all considered 
3d, 4d (not shown) and 5d TMs (Fig. 2(c)). For example, 
the break forces with and without spin-polarization for 
W differ by a factor of two. Both at the beginning and 
the end of the 5d series, this ratio reduces according to 
smaller maximal magnetic moments and is equal to one 
for the non-magnetic noble metal Au. The break force 
Fq not only represents the maximal force applicable to 
a chain, but also serves as a measure for the maximal 
amount of mechanical energy which can be gained upon 
relaxation of its bonds. In any case, high values of Fq 
favor successful chain creation in BJ experiments. Thus, 
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FIG. 2: (color online) (a) Cohesion energy difference A.Ehead, 
(b) difference d — do = In 2/7 and (c) break force Fo for 
non-magnetic (marked as NM, black squares) and magnetic 
(marked as FM or AFM, green circles) 5d-TM chains. Broken 
line in (c) stands for the break force Fq 1 calculated accord- 
ing to Eq. ©. As indicated by the grey shaded area, W 
and Re MWs reveal an AFM ground state at all inter-atomic 
distances, while the rest of the chains are FM. With orange 
squares in (c) the break force for the FM W and Re MWs is 
shown for comparison. 



as the appearance of magnetism crucially reduces Fq, we 
can conclude, that the formation of local magnetic mo- 
ments suppresses chain formation. 

In order to capture the origin of the magnetically in- 
duced reduction of Frj, we relate the magnetic quantities 
to the non-magnetic ones. In the first approximation, the 
influence of the magnetization on the binding energy po- 
tential £ sp (d) = £nm {d) ~ £m (d) can be attributed to the 
Hund-type intra-atomic exchange interaction ^IM 2 (d) 
between mostly d- but also s-electrons on the atomic 
sites and the inter-atomic Heisenberg J(d)Mi(d)M i+ i(d) 
exchange contributions between the atomic spins. The 
binding energy £m(oo), prefactor in Eq. (1), then reads 
£m(oo) = £nm(oo) - fsp(oo) + £ sp {d ), where £ sp (oo) = 
i/A/ 2 (oo) is purely given by intra-atomic exchange of 
a free atom and £ Bp (do) = £nm(c?q M ) — £M(d^) is the 
energy difference of the NM and M states at the corre- 
sponding equilibrium distance. Taking as an example 5c? 
TMs we can safely assume even for W, exhibiting the 
largest equilibrium magnetic moment of M(d^f) « 1/ib 
through the series, that the impact of magnetism on 
the equilibrium properties of the chains is small and set 



for simplicity d^ 1 = d$ M and £ sp (do) = 0. As appar- 
ent from Fig. 2(b), it is also reasonable to assume that 
7m = Tnmj in which case the magnetic break force sim- 
plifies to: 



* = — -£m{oo) = Fq — -I d M (00). (2) 



For 5d TM chains we estimated F^ according to Eq. ([2]) 
using ab initio values of the non-magnetic break force 
Fq M , magnetic moments at d = 6.5 a.u. and the atomic 
exchange integrals between ri-electrons, Id, from Ref. [27I ] . 
and plotted it in comparison to the magnetic break force 
determined from ab initio in Fig. 2(c). Good qualitative 
agreement between the two break force curves underlines 
that the intra-atomic exchange / is the major origin of 
the magnetically driven reduction of Fq, while the inter- 
atomic exchange plays only a minor role. This conclusion 
is further verified by the observation that changing the 
magnetic order from AFM to FM in W and Re MWs 
results only in a small change of F (see Fig. Efc)). 

We now turn back to the analysis of our ab initio re- 
sults. With the knowledge of all key quantities entering 
the criteria for stability and producibility we can further 
analyze both criteria in the phase space of the number of 
atoms N and inter-atomic distance d presented in Fig. 3. 
Each of the criteria leads to a distinct region where it 
is fulfilled and accordingly the chain is stable (S) or pro- 
ducible (P). Ideally for a successful chain elongation event 
to happen, both regions (S) and (P) have to overlap (SP). 

Comparison of the phase diagrams for different TMs, 
shown for Fe, Ru, and W in Fig. 3, underlines, that the 
formation of local magnetic moments strongly suppresses 
the probability of chain formation for 3c?, 4c?, and 5c? el- 
ements. If we ignore the formation of magnetism among 
the 3c?-, 4c?-, and 5c?-TM series Cr, Mn, Fc, Ru, Rh, 
Ag, Re, Os, Ir, Pt and Au exhibit extensive SP regions 
(Fig. 3(b), (d)) indicating successful chain formation for 
these elements. Even for W (Fig. 3(f)) with bcc(110) 
electrodes neglecting magnetism results in touching S and 
P regions, indicating chain formation for more open lead 
structures. Allowing for the formation of local spin mo- 
ments the picture changes completely: SP regions emerge 
exclusively for Ru (Fig. 3(c)), Ag, Ir, Pt, and Au while 
for all other elements the S and P regions are clearly 
separated and no chain formation occurs (Fig. 3(a),(e) 
and in [3]). While for Pd and Pt chains the influence of 
magnetism is small due to relatively small moments en- 
tering Eq. ([2|), the SP-regions of Ru (Fig. 3(c)) and Ir arc 
considerably less extended for magnetic chains than for 
non-magnetic ones, underlining the suppression of chain 
formation by magnetism. 

While our predictions based on the assumption that 
chains in BJs are magnetic match and explain the exper- 
imental findings for successful Ag, Ir, Pt and Au chain 
formation, the results of the model for non-magnetic sus- 
pended chains contradict the experimental observations 
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FIG. 3: (color online) Phase diagrams for Fe, Ru, and W- 
BJs with (a),(c),(e) and without (b),(d),(f) spin-polarization 
of the chain atoms. Plots indicate regions of stability (S, dark 
gray), producibility (P, light gray), separated by a white re- 
gion or overlapping in the SP region (in green). The starting 
point along the x-axis do and the assumed surface orienta- 
tions of the leads are indicated. The input parameters for 
W are given in Fig. 2(a)-(c) and for Fe and Ru in [2(|. (g) 
provides a schematic summary of our results. Shown is the 
chain formation probability, proportional to the size of SP- 
regions, versus chain length, proportional to the highest N 
for which an SP-region exists, with (green) and without spin- 
polarization (gray) in relation to the experimental findings 
(red) [ToL [l3l . [l4 . fig . Il8l . [29| in all cases in arbitrary units rel- 
ative to Au. Arrows indicate (exemplified for Fe, Ru and W) 
the consequence on the chain formation caused by switching 
off the finite magnetization in the chain. 



at several crucial points (Fig. 3(g)). Firstly, nanocontacts 
of 3d-TMs such as Fe are reported to form only point- 
contacts with no tendency to form longer chains [l6| . 
Secondly, also BJ-experiments using W as tip-material 
result only in point-contacts [28j|, moreover, W tips are 
widely used in STM- and AFM-experiments due to their 
structural rigidity preventing substrate-induced reforma- 
tions [29l | . Thirdly, non-magnetic Ir chains would become 
as long as those of Pt, and almost as long as those of 
Au, in direct disagreement with experiments which re- 
port significant decrease of chain formation probability 



and length when going from Au to Ir [lfj, |13|, [LJ, 1 1 81 ] - 

These clear contradictions to existing experimental evi- 
dence lead us to the conclusion that only when chains 
in BJs are magnetic, the experimentally observed trends 
can be reproduced and explained throughout the 3d, Ad, 
and 5d transition- and noble-metal series. Therefore, by 
reductio ad contradictum, comparing theoretical predic- 
tions with experimental findings, we provide a convincing 
evidence that TM chains in BJs are magnetic. 

We thank J. M. van Ruitenbeek for fruitful discus- 
sions. S.H. thanks the Stiftcrvcrband fur die Deutsche 
Wissenschaft for financial support. 
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